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Abstract 

We investigate the quantum evolution of the metric operators for Bianchi-Type I model 
universes in the Heisenberg picture in order to remove the need to consider the wave function 
of the universe and interpret its "spin" variables. The calculation is analogous to that of 
the Zitterbewegung of the Dirac electron. We consider the behavior of the metric near the 
classical singularity, and consider the curvature there. Although factor ordering questions 
preclude the presentation of an unambiguous result for the curvature invariants, it does seem 
that the classical t~ 4 divergence of the Kretschmann scalar is not removed by quantization. 



1 rbmann@sciborg.uwaterloo.ca 
2 email: gmurphy@raex.com 



1 Introduction 



The Hamiltonian treatment of simple cosmological models has provided one useful approach to 
quantum cosmology. In the classical case the canonical equations have solutions which represent 
solutions of Einstein's field equations, and the hope is that the corresponding quantum equations 
will represent a full and consistent theory of quantized gravitation [1]. The Hamiltonians for 
homogeneous universes in general relativity occur in squared form, the simplest example being 
that for empty Bianchi-Type I spaces with two degrees of freedom: 

h 2 = p 2 + +p 2 . (1) 

where we take c, but not h to be one. 

The problem of quantizing a system with such a Hamiltonian already had to be faced in the 
early days of relativistic quantum theory. The answer which led to the greatest progress was 
that of Dirac, who extracted the square-root of the corresponding expression for a particle by 
introducing spinors. With H 2 = p 2 + m 2 , one writes 

H — a • p + (5m (2) 

with (a,/3) a set of anti-commuting 4x4 matrices. This led, as is well known, to a natural 
understanding of electron spin and related phenomena. 

The same type of procedure can be carried out for cosmological models with a few degrees of 
freedom. With the Hamiltonian given by (1) we write 

H = a + p+ + (3) 

where the cr's are 2x2 Pauli Matrices. Plane wave solutions of the resulting Schrodinger equations 
are then easily found. But this procedure, while mathematically quite simple, raises interpretative 
problems. The most obvious difficulty is that the significance of states corresponding to "spin up" 
and "spin down" for a cosmological model remains obscure [2]. 

Since this difficulty has to do with interpretation of the state vector, it may be helpful to 
follow a route that focuses instead on the dynamical variables of the model. Rather than work in 
the Schrodinger picture, as has usually been done in quantum cosmology, we may deal with the 
problem in the Heisenberg picture. This approach provided some novel insights into the dynamics 
of the Dirac electron, and may be expected to do the same for a model universe. 

2 The Zitterbewegung of the Electron 

For the Dirac Hamiltonian (2), the Heisenberg equations of motion give the velocity of the particle 
as dx l jdt = [x\ H]/ih = a 1 . Proceeding further, we obtain ifida 1 jdt = —2a l H — 2p % , a differential 
equation with the solution a 1 = A % exp(—2iHt/H)+p l H~ 1 , where the Ai are constant. Since a 1 = 1 
and p l H~ l < 1, the magnitude of the A 1 must be of order unity. Another integration yields 

x i = {ihA 1 /^) exp(-2iHt/h)H- 2 + ^H'H + a 1 (4) 
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with the a 1 further constants. The second and third terms in (4) represent the motion of a particle 
that is expected from classical dynamics, motion in a straight line with constant velocity p/H. 
The first term, on the other hand, represents a quantum oscillation about that classical motion 
with high frequency (> 2m/ H) and an amplitude comparable with the particle's reduced Compton 
wavelength. This is the Zitterbewegung[3] or "trembling motion". When the Dirac equation is 
taken to describe a single particle, the peculiarly quantum motion provides one way of thinking 
about the spin of that particle. 

3 Quantum Cosmology in the Heisenberg Picture: The 
Metric 

In the cosmological problem, the co-ordinates canonical to p are denoted by f3, and we use a time 
parameter f2, which is related to the conventional cosmic time t by t — (4irR3/H)e~ 3n , where R 
is a constant. (Thus t — > corresponds to VI — > oo.) The diagonal spatial metric is given by 

gn = R 2 e- 2n exp[2(P + + V3(3-)\ 

g 22 = R 2 e- 2n expW+-V3P-)} (5) 
g 33 = i? 2 e- 2n exp[-4/3 + ] 

The classical equations of motion with (1) then give simply (3± = (p±/H)Q or 

/3+ = Qcos0, /3_ = ftsm0 (6) 

if /3(0) = 0, where cos# = p + /H and sin# = p^/H. When these expressions are substituted into 
(5), we obtain a representation of the well-known Kasner metric[4]. 

We now move to the corresponding quantum calculation with the Hamiltonian (3). If we use 
a + a_ + a + a_ = 0, o"^ = o" 2 = 1, the Heisenberg equations of motion give us, in analogy with (4), 

(3± = (ihB±/2) exp(-2iHn/h)H- 2 + p ± H~ l Vt + C± (7) 

Here the S's and C's are constant 2x2 matrices. 

It is tempting to think that the C's can essentially be ignored. A similar constant of integration 
representing the initial position of a particle can often be dropped in elementary mechanics by just 
defining the initial position to be at the origin of coordinates. But in quantum theory matters are 
different, as we can see by considering the elementary problem of the motion of a non-relativistic 
free particle with H = p 2 /2m. 

The Heisenberg equations of motion easily give p = constant and q = (p/m)t + go • We might 
be tempted to ignore the constant go but it cannot vanish because the commutation relations 
between g and p require that [go,p] = ifi- 

Having to retain the C's as additive constants might not seem to be a serious problem. But in 
our case the required exponentiation of the /3's in the metric mean that we will have multiplicative 
constant operators which do not commute with the p's in the metric. This will complicate the 
quantum calculations. 



2 



The metric components, obtained by substituting (7) into (5), are now operators which carry 
the time-dependence of the model. They are 

0ii = {expDi)R 2 exp ih{B + + V3B_) exp(-2iHVt/h)H- 1 exp 2(cos6 + V?, sin# — 1)0 
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exp D 2 )R exp \ih(B+ - V?>B^)exp(-2iHVt/K)H- l \ exp 2 (cos 9- v^sinfl-ljO 



g 33 = (exp D 3 )R 2 exp [-2iHB + exp(-2iHVl/K)H~ l ] exp [-20(1 + 2 cos 0)] 

where the exponentiated -D's are the multiplicative constant operators noted above. 
In terms of the cosmic time t we have 



0ii = (exp D,)R 2 exp\iH(B + + V3B^(Ht/S) 2iH / 3H H- 1 ] (/ft/s)2(cos0W3sin0-i)/3 

22 = (expD 2 )R 2 exp \ih(B+- y/ZB_){Ht/ Sf^'^H' 1 ] ( H t/S) 2{cosd -^ sin9 - 1)/? > (9) 
033 = (expD 3 )R 2 exp[2thB4Ht/S) 2lH/m H~ 1 ](Ht/S) {1+2cosd) 

where S = AttR 3 . The third exponential factors in (8) and (9) give the classical Kasner behavior, 
while the preceding factors contain the distinctive quantum mechanical behavior. It should be 
noted, though, that cos# and sin#, now defined by (6) as p + H^ 1 and p_H~ l respectively, are 
operators, so that the expressions for do not split cleanly into "classical" and "quantum" 
factors. 

Sinceif 2 = p 2 = p\ + p 2 __ we have exp(aH) = cosh(ap) + (H/p) sinh(ap). If we use this result 
in (7) we find 

(3 ± = [(H/2p 2 )B ± sm^)+p ± Q/p](H/p) + ( i h/2p 2 )B ± cos^) + C ± 

The metric components can then be written as 

gu = (expD i )R 2 exp[\ i (t) +'y i (t)H/p] 

where 



(10) 



A 3 

7± 

73 



-20 + i[H(B+ ± VSB_)/p 2 ] cos(^) 



-20 - 2i(hB + /p 2 ) cos' 



20p, 



= (h/p 2 )(B + ± V3B_) sin(^) + 2{p + ± v / 3p-)0/p 
= -2[(H/p 2 )B + sm(^)+2p + n/p] 



(11) 



with D 1 = 2(C+ + v/3C_) , D 2 = 2(C + - v/3C_) , and D 3 = -4C+ . (We also have dropped the 
summation convention.) 
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For brevity we put fa = \ + ^iH/p so that ga = exp (Dj) R 2 exp (<pi). The contravariant 
components g w = R~ 2 exp (—fa) exp (— Dj) will then be the correct inverse of ga . 

The temporal evolution of the metric is now rather complicated. The quantum feature par 
excellence is the fact that B and H are operators which can be represented by 2 x 2 matrices 
in the present model. But we should also note that the simple complex exponential or circular 
functions which describe the Zitterbewegung of a particle are now themselves the arguments of the 
exponential functions which give the metric in (8) and (9). As t — > and Q — > oo, the magnitude 
of the metric components will oscillate an infinite number of times. 



4 Quantum Cosmology in the Heisenberg Picture: The 
Curvature 

In order to investigate the behavior of space-time near t — (Q — oo), the true singularity of the 
classical model, we must of course examine not simply the metric but the curvature. 

In the classical case we can use the cosmic time t and write the metric corresponding to (9) 
as <?oo — — 1) 9aa = t 2sa with appropriate scaling. The no n- vanishing components of the Riemann 
tensor are given by R a0a0 = -r 0a a,o + r aa0 r 0aa = -s a (s a - l)t 2sa ~ 2 and the Kretschmann scalar 

K = _16^i£H£i (12) 
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which diverges as t — > 0. 

The quantum theory is more complicated and there is considerable ambiguity because of the 
non-commutativity of operators. There are three ways in which this complicates matters. First, 
there is the factor-ordering problem, which always arises when a transition is made from a classical 
to a quantum theory. Secondly, the operator scale factors in the metric which we noted above 
must be taken into account. Finally, we have to be aware that the B's must have an algebraic 
structure ensuring that cr's (which now in the Heisenberg picture are operators) have the correct 
algebra. 

If we choose the usual order of factors for the Christoffel symbols we obtain 

X i 

r m0 „ = ^mn^'m, Omn = --Ro( eX P D m )5 mn (f)' m (exp <f> m ), 

where a prime denotes the time derivative. We then calculate curvature components with the 
ordering of factors given by 

E)A _ pA pA _|_ per pA per pA /-i o\ 

)J,UK (IV, K UK,U ' flU K<J UK VG \ °) 

It would be desirable, as we will see below, for factors to be ordered in such a way that all 
occurrences of exp D m in the curvature occur to the left of other factors. Unfortunately this 
would require a different ordering of factors for different components, so that such a result would 
not be covariant. Instead of trying to indicate all the possibilities we will simply follow the 
prescription in (13) and write the results as follows. 
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R° 



R l 



mnO 



rank 



exp D m (Rl)5 n 



<t>m + 



exp0 m + C.T 



(14) 



expD m (Rl) {5 mn 5 k i - 5 mk 5 n i) 0' m 0; + C.T 



Here C.T. denotes "commutator terms," distinctively quantum mechanical expressions involving 
fi, which result from switching the order of factors containing D's (and thus C's) and p's. 
We then find that the mixed Ricci tensor components are 



R 



3p 2 f 



B 2 



cos 



2VLp 



+2pBcos ( ) 2/ sin 



Rl = 



Rl 



Rl 



4i 




-2i 



cos 



COS 



— sin 

2Qp 

~h 

2Qp 

~h 

2Qp 



2VLp 



(^1 • pB 



+ i sin 



+ i sin 



h 

2ttp\ H 
h J p 

2ttp\ H 



(J^) [ cos ( 



IT 



+ i sin 



« y v n 

2VLp\ H 
hip 



P 



+ C.T 
+ C.T 



(15) 



+ C.T 



Because of cancellations the Ricci scalar R is simply equal to Rq. We have written pB =p + B + + 
p_B_ and B 2 = B\ J r B 2 _. In the calculation of the mixed Ricci components, all occurrences of 
the commutator of B + and B_ cancel out. 

Calculation of the Kretschmann scalar is now considerably more complicated, even if the 
problems of factor ordering are ignored. In addition, the need for four factors of metric components 
to lower and raise indices introduces a great many more choices for their ordering, and the non- 
commutativity of the Bs can no longer be ignored. The full expression for K with the choices we 
have made for factor ordering is quite lengthy. It can be written as 



H 

K = K + K 1 — 
p 



(16) 



where 



Kn = 



27hH 4 



IB 2 ( 1 - 2 cos 2 ( ^p ) ) + 5h 2 



+ 



32 sin( 2 fp) 



27Ht 4 p 2 



2m A p 

[-3H (B 2 - p 2 ) + 8i (B 2 +P+ - B 2 _p + - 2p_B + B_)] 



, 2Q . 
25^ cos ( — p ) + p-B 
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^i 6 ! „ cos (^P) {Sihp-B-QB + (b+ + V3B^ (b+ - y/ZB^j + 2 (p 2 + B + - p 2 _B + - 2B_p + pJ) } 



16* cos 2 (fp) 



27Ht 4 p 2 
128i cos 3 (^p) 



27fa 4 p 2 



[3i/LB 2 + 8 (5 2 p+ - - 2p_ 

-£+ (5 + + v^S-) (b + - v 7 ^-) 



12g sin(fp) cos 2 (fp) 

27tV 
64isin(fp) /2Q 



B + [B + + VSB-) [B+- V3B 



9tV 
32i sin(fp) 



(-6fi_p + p_ + 3B+B 2 + 3B + p 2 + - 3B + p 2 _ - Bl) 

cos 3 (j^-p) (5 3 p- + B + B 2 _p + + 5 3 p + + B 2 + B_pJ) 



27tV 

32 A fm \ _64_ 



+ 3^4 c ^ 2 (^V) (5*p 2 _ + 35 2 p 2 - 25 4 + 35 2 p 2 + AB_p^B + p + - AB 2 + B 2 _ + B 2 _p 2 + ) 

+ ^7 cos (7^) + Pl B + + P+P-B+ ~ 3B\B_p_ - 3B%p + - 3B 3 _ P - - 3B + B 2 p + + B_p\p_) 



+ ^i(p 4 + BA ~ 2B lP- ~ " QB 2 p 2 _ - 8B_p_B +P+ - 2B 2 _p\) 



(17) 



64i 



2Q 



27tw sm {T p ) C0S {T p ) B + 27m P 



2Q 



32i 



B< 



, , 2Q 
1—2 cos - I — p 



'20, 

2p ■ B cos ( — p 



128 sin (^p) cos 2 (^p) / / /- \ / /- \\ 

16sin (f p) cos m {3ifiB * + 8 ( s 2 p + - 5 2 p + - 2p_3 + 3_)) 



27fa 4 p 2 
16sin(fp) 



27Ht 4 p 2 



128 cos 3 (fp 



ft - 



27*V 
32cos(fp) 



(3ifip • 5 + 2 (p + 5+ - p 2 fi + - 2p_p+B_) + 65+fi 2 - 2fi 3 ) 

[B + (B + + V3S_) - V3S_)) - 64c ^(f p ) _ S 2 p+ _ 2p _ s+s _) 



9tV 



(3S+5 2 - 5 3 + (p 2 + 5 + - p 2 B + - 2p_p + 5_)) 



32 



27t 4 p : 



■{QB + B_p_ + 3B 2 p+ + p 3 - 3B 2 +P+ - 3p 2 _p+) 



(18) 



32isin (fp) cos 3 (fp) a 64i sin (fp) cos 2 (fp) , „ , , 9 

+ V ^ ; ^-iB 4 + V * *7 K -^{p + B + B 2 _ + + 5 3 p + + 3B 2 + B_p_) 



3t 4 p 2 



3tV 
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16i sin (^p) cos (^p) 
3tV 



(-2B 2 + B 2 _ + 4B+£_p_p+ + 3p 2 5 2 + 3^5 2 + 5 2 p 2 - fi 4 + fi 2 ^ - 5+) 




(- P+ Bl + p 2 + p 3 + 5 + + p\p-B + 



-B 



B 2 _p+ - B 2 + B_p_ + B_p 3 _ - BlpJ) 



The terms K that do not involve any distinctively quantum factors, and that are not affected by 
factor ordering complications, can be written as 



This is identical with the classical result given by (12). This means that while the magnitude of 
this curvature invariant will undergo oscillations of increasing rapidity as t — > 0, there will be an 
overall unbounded increase beyond any bounds according to the same t~ A law that obtains in the 
classical case. The leading divergence of the curvature is not removed by quantization. 

There are, of course, profound differences between classical and quantum curvatures in other 
regards. As is well known, the special case of the Kasner metric with s 1 = s 2 = and s 3 = 1 in 
the classical theory is simply flat space-time. This is no longer so in quantum theory because of 
the additional exponential factors in the expressions (9) for the metric. 

Use of the Heisenberg picture thus gives some novel insight into the quantum dynamics of the 
simplest model universe of the Bianchi type, and shows in particular that the initial singularity 
is not avoided by quantization. There are no questions about boundary conditions on a wave 
function or the interpretation of it to cloud this result. Unfortunately the procedure followed here 
cannot be extended immediately to more complicated models. Just as straightforward integration 
of the equations of motion for the Dirac electron works only for a free particle and not when 
an external potential is present, the curvature term which would be added to right side of (1) 
frustrates attempts to follow the same procedure for spaces of other Bianchi types. 

In a sense our result is not too surprising. When we solve the Dirac eqn for a particle we 
recover the classical motion of the particle together with the distinctive Zitterbewegung effects. 
In the cosmological case a similar thing happens - we recover the classical cosmological motion 
corrected by quantum effects. But in this case the classical solution is divergent, and the question 
of interest is whether or not the quantum effects can cancel this divergence. In particular, does 
the expectation value of the Kretschmann scalar blow up as t — > in the quantum case as it does 
in the classical one? 

We have not been able to answer this question definitively because of the huge number of 
possible factor orderings in K that arise in different ways. The most straightforward way of 
writing K does give a divergent result, and this may suggest that hopes for an elimination of the 
initial cosmological singularity by quantum theory will be disappointed. But it is not impossible 
that some choice of factor ordering will eliminate this divergence. 

Use of the Heisenberg picture thus gives some novel insight into the quantum dynamics of the 
simplest model universe of the Bianchi type, and shows in particular that the initial singularity 
is not avoided by quantization. There are no questions about boundary conditions on a wave 
function or the interpretation of it to cloud this result. Unfortunately the procedure followed here 
cannot be extended immediately to more complicated models. Just as straightforward integration 




(19) 
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of the equations of motion for the Dirac electron works only for a free particle and not when an 
external potential is present, a term involving the 3-space curvature term which would be added 
to right side of (1) frustrates attempts to follow the same procedure for spaces of other Bianchi 
types. The problem would, of course, be even more difficult for nonhomogeneous universes or 
other space-times [5]. 
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